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^! Abstract 

The concept of quantum speed limit-time (QSL) was initially introduced as a lower bound 
■ to the time interval that a given initial state ipi may need so as to evolve into a state or- 

^ . thogonal to itself. Recently [V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. A 67, 

\ 052109 (2003)] this bound has been generalized to the case where "0/ does not necessarily 

' evolve into an orthogonal state, but into any other ^p. It was pointed out that, for certain 

o ■ 

\^ • classes of states, quantum entanglement enhances the evolution "speed" of composite quan- 



o 



tum systems. In this work we provide an exhaustive and systematic QSL study for pure 
and mixed states belonging to the whole 15-dimensional space of two qubits, with ipp a not 



■ necessarily orthogonal state to ■0/- We display convincing evidence for a clear correlation 

^3 : 

O^. between concurrence, on the one hand, and the speed of quantum evolution determined by 

. ^ the action of a rather general local Hamiltonian, on the other one. 
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I. INTRODUCTION 



One of the most fundamental concepts in the quantum description of nature is 
that of entanglement , which in recent years has been the subject of 



intense research 
references: 




the following, by no means exhaustive list of 
A state of a composite quantum system is 
called "entangled" if it cannot be represented as a mixture of factorizable pure states. 
Otherwise, the state is called separable. Entanglement constitutes a physical resource 
that lies at the heart of important information processes such as quantum 

teleportation, superdense coding, and quantum computation. 

Entanglement is essential for both (i) our basic understanding of quantum mechanics 
and (ii) some of its most revolutionary (possible) technical applications. Thus, it is 
imperative to investigate in detail the relationships between entanglement and other 
aspects of quantum mechanics. In such a vein, particular interest is assigned to the 
exploration of the role played by entanglement with regards to the dynamical evolution 
of composite quantum systems. 

In this effort we will be interested in the speed up of quantum evolution produced 
by entanglement. Why? Because in quantum computation one tries to (i) avoid loss 
of coherence and (ii) increase the velocity of information processing and information 



transmission. Anandan and Aharonov 



ll\ have shown that, given a state \ip) and a 



curve C in the projective Hilbert space V, the quantity 



dt--, (1) 



with 



{AEf = (m'm - imm', (2) 

is independent of the particular Hamiltonian H used to transport the state along the 
curve and is in fact the distance along C as measured by the Fubini-Study metric, 
deducing as a consequence the uncertainty relation 

{AE)At > ^, (3) 

where (AE) is the time-averaged uncertainty in energy during the time interval At. 
Equality in Eq. ^ holds iff the system moves along a geodesic in V. In this case the 
evolution may be said to have minimum uncertainty, analogous to how a Gaussian wave 
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packet is said to have minimum position-momentum uncertainty at a given time. More 
generally they define an efficiency in evolution e = s/sq, where the denominator gives 
the distance along the shortest geodesic joining the initial and final points of evolution. 
Loss of coherence in evolution may be regarded as due to the time-energy uncertainty 
principle. In trying to avoid such loss, speeding up evolution seem advisable. 

In this regard, Margolus and Levitin have shown that the minimum evolution 



time in which one state evolves to an orthogonal one depends on the mean energy 
and the fiuctuation. Giovannetti, Lloyd, and Maccone 0] recently uncover the 
fact that, in certain cases, entanglement helps to "speed up" the time evolution of 
composite systems. This "speed" of quantum evolution is also of considerable interest 
because of its relevance in connection with the physical limits imposed by the basic 
laws of quantum mechanics on the velocity of information processing and information 
transmission [3, ll^ ll^l. 

The evolution "speeding-up" ability of entanglement has been demonstrated only 
in special, if important, instances. One would like to ascertain that it is indeed an 
entanglement feature, and not just something that happens in these instances. Thus 
we will here undertake a general study. The aim of the present contribution is to make 
a systematic study of the connection between (a) entanglement and (b) the speed of 
quantum evolution as determined by the action of a rather general local Hamiltonian, 
by means of a numerical simulation. Our model belongs to a family that includes the 
)asic models of quantum optics and cavity QED [n a previous work 

2^ a corresponding study was performed just for (i) pure states of bipartite systems 
of low dimensionality evolving towards (ii) an orthogonal state. Two different cases 
were analyzed: i) two qubits (distinguishable) systems and ii) bosonic or fermionic 
composite (bipartite) systems of the lowest dimensionality. In the present effort we 
are going to tackle an extension to the case of two distinguishable systems of (a) pure 
and (b) mixed states that evolve to (c) any other state, not necessarily orthogonal to 
the initial one. We also consider (d) the special case of maximum entangled mixed 

n 

states (MEMS) [21], and also (e) that of the set of mixed states whose entanglement 
degree cannot be increased by the action of quantum gates (IH states) |22| . We remark 
on the facts that i) MEMS have recently been detected experiment allyj^, and ii) 
nowadays the possibility of obtai ning such states via the action of local non-unitary 



quantum channels is being studied 



25 



26|. Thus, the ensuing results will be applicable 



to any physical systems where bipartite states play a leading role. 

The paper is organized as follows: in Sec. II, using the time evolution of the 
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fidelity we present the quantum speed limit for pure states. The case of mixed states 
is presented in Sec. III. The special case of MEMS and IH states are presented in Sec. 
IV, and finally some conclusions are drawn in Sec. V. 

II. QUANTUM SPEED LIMIT FOR PURE STATES 

Let us consider first the dynamical evolution of pure states for the case of two equal 
but distinguishable subsystems evolving under a local Hamiltonian, that is, we deal 
with a two-qubits system whose evolution is governed by the (local) Hamiltonian 

H = Ha®Ib + Ia<^ Hb, (4) 
whose eigenvalue equation writes 



Ha,b = £a,b\^), 
Ha,b = 0|0), 
= ^ac, 

eB = Sbg, (5) 
e being an arbitrary energy. 

Our bipartite states (the eigenstates of H) are |00), |01), |10), and |11), while the 
concomitant eigenvalues equal 0, 6b€, 6ag, and {6a + ^b)^, respectively. 

A. General methodological considerations 

In this paper we perform a systematic numerical survey of the evolution properties 
of arbitrary (pure and mixed) states of a two-qubits quantum system, under the action 
of the Hamiltonian by recourse to an exhaustive exploration of the concomitant 
state-space S. To such an end it is necessary to introduce an appropriate measure fi 
on this space. Such a measure is needed to compute volumes within S, as well as to 
determine what is to be understood by a uniform distribution of states on S. The 
measure that we are going to adopt here is taken from the work of Zyczkowski et al. 
js^, I3| • An arbitrary (pure or mixed) state p of a quantum system described by an 
A^-dimensional Hilbert space can always be expressed as the product of three matrices, 

p = UD[{\}]U\ (6) 
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Here ?7 is an x unitary matrix and i5[{Aj}] is an x diagonal matrix whose 
diagonal elements are {Ai, . . . , Xn}, with < Aj < 1, and J2i -^j = 1- The group of 
unitary matrices U{N) is endowed with a unique, uniform measure: the Haar mea- 
sure u j^l- On the other hand, the A^-simplex A, consisting of all the real A^-uples 
{Ai, . . . , Aat} appearing in Eq. ©, is a subset of a (A^ — l)-dimensional hyperplane 
of 72.^. Consequently, the standard normalized Lebesgue measure Cn-i on IZ^^^ pro- 
vides a measure for A. The aforementioned measures on U{N) and A lead then to a 
measure fi on the set S of all the states of our quantum system 21, 2^,12^, namely. 



= vCn-i- (7) 

In our numerical computations we randomly generate pure and mixed states ac- 
cording to the measure (jZj). 



B. Pure states 



For pure states |\1') of our composite system the natural measure of entanglement 
is the usual reduced von Neumann entropy S[pa,b\ = —TrA,B{pA,B^^PA,B) (of ei- 
ther particle A or particle B) where pa,b = TrB,A{\'^){'^\)- It is convenient for our 
present purposes to use, instead of information measure S[pa,b] itself, the closely re- 
lated concurrence value C, given by 

C^ = 4detpA,B. (8) 

Both the entanglement entropy S[pa,b] and the concurrence C are preserved under the 
time evolution determined by the local Hamiltonian (^. Given an initial state 

\^it = 0)) = ColOO) + CilOl) + C2IIO) + C3III), (9) 

with 

4 

Y,\ci\' = l, (10) 

2=1 

its concurrence is. 



= 4|coC3-CiC2|l (11) 

Our objective is to characterize the departure of the system, at a time t (represented 
by \E'(t)), from its initial state \E'(t = 0). To this end we can use the quantum concept 
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of fidelity P that, for pure states, is the squared- modulus of the overlap between the 
two states involved, i.e.. 



P{z) = \{^{t = 0)|^(t))| 2 = I |co|' + |ci|V^ + |C2|V-^ + |c3|V-^+^^ I 2, (12) 
where 

z = exp(ifi), (13) 

and 

a^l (14) 

From now on we consider time intervals measured in units of h/e, and use for this 
rescaled time the letter Q. We will also use the following notation: = 6a + Sb and 

A- = 6a- Sb 

The key idea is that of measuring the speed of dynamical evolution by studying the 
time evolution of the fidelity. To such an end one first of all fixes a given P amount, 
say P = F, and proceeds to calculate amount of time needed for a given state to evolve 
from P{z) = 1 (at t = 0) to P{z) = F at, say, t = r, for F G [0, 1]. 

In Ref. only the case 6a = 5b = ^ was discussed. This particular Hamiltonian 
instance will be referred to as corresponding to the Hamiltonian Hj in Sec. IV. 
The condition (fT^ speciahzes for Hj to 

F = 2po3cos2fi + 2po3(l-so3)cos^] + (l-so3)' + |co|' + |c3|^ (15) 

where pos = |conc3p and so3 = |coP + |c3p. 

In this case a minimum of the fidelity is achieved for the special value Q = Qmin 
given by 

n ~ '503)S03 

ilmin = arccos , (16) 

4po3 

which can yield nonphysical complex values. To avoid this we limit the argument of 
the arccos to the interval [—1,0]. 

The time r required to evolve to a state with fidelity F admits a lower bound that 
depends upon both the state's expectation energy E and its fluctuation AE ||^, 

TL.Bound = maxia{F) — ,p{F)^^y (17) 
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where the functions a{F) and P{F) are detailed in Ref. 121. We can compute P{F) 
using an expression previously proved in Refs. 0| 



I3{F) = - arccos 



while a{F) can be numerically calculated with great accuracy For F = (when 
the initial state evolves to an orthogonal one), one finds <y{F) = f3{F) = 1. In the 
opposite situation, when the state does not evolve (F = 1), both functions vanish. 

n 

In Ref. j20|| situations were dealt with for which the orthogonal state to the initial 
one was definitely reached, which is not the usual case. A useful parametrization, 
introduced in this reference, reads 



I |2 I |2 -p 

|co| = icsi = r, 

|ci|^ = — 25rcosi7 

\c2\^ = -2(1 - (5)rcosfi, (19) 

with r = 1/[2(1 - cos^])] and n E [f,7r], with S E [0,1]. In other words, = 
arccos((2r — l)/2r). We note that introduction of this parametrization in the right- 
hand-side of Eq. (fT^ for F = (orthogonality) does yield an equality. 

To study the role of the entanglement on the speed of quantum evolution and on its 
lower bounds, one should pay attention to the {C,t /Tl. Bound) plane. A representative 
group of two-qubit states evolving to an orthogonal one is depicted in Fig. 1(a). These 
states can easily be generated using the parametrization (fT^ . The ratio t/T^ Bound has 
a maximum value equal to no matter which value C adopts. The minimum value 
of this ratio does strongly depend on C through the value of Q, and can be analytically 
obtained [2^. Only maximally entangled states reach the bound Tl. Bound- Separable 
states have a different behavior. For a rather general Hamiltonian of type Hj, they all 
evolve to an orthogonal state in a fixed time t/Tl Bound = v^- 

These features can be easily explained. For pure states evolving to an orthogonal 
state according to if/, the minimum time interval required to complete such evolution 
depends only on Ai?. This is the first option in Eq. (|T7j) . In this specific case AE 
strictly depends on the value of Vt [2^ . All separable states evolve to an orthogonal one 
in a fixed time = tt, for which the minimum possible value of the ratio t/Tl, Bound is 
\/2. Thus, for separable states the maximum and the minimum of such ratio coincide. 
As the entanglement is increased, pure states are able to evolve in more rapid fashion. 
This is the reason that lies behind the dependence of the minimum of t /T^ Bound with 
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FIG. 1: t/Tl. Bound for pure states, a) Pure states that evolve to an orthogonal one and b) 
pure states for which F^m G [0.35,0.4]. The points Pi and P2 represent the fastest states 
corresponding to each of the families determined by Eq. (|17|) . For the range of values of 
Fmin considered here, these fast states correspond to Fmm = 0.35. 

the concurrence. Maximally entangled states can evolve to an orthogonal state in any 
time lapse within the range Vt G (^, tt). The states that evolve in the shortest possible 
time {Vt = vr/2) are also those to which a minimum of the ratio t/T^ Bound = 1 is 
assigned. Of course, we cannot extend the same conclusions to the rest of the states 
that evolve to some value of F 7^ 0, as it is clearly seen in Fig. 1(b), or to other 
Hamiltonians than Hj. 

From Eq. (fT^ we realize that the fidelity for pure states oscillates in time. Thus, as 
time goes on the fidelity of a given state climbs and goes down in alternating fashion, 
reaching minima of different depths. Our interest in this respect will be focused on 
those special times at which the first fidelity-minimum Fmin is attained. Doing so we 
can treat all pure states in a unified manner. Note that the type of oscillation we are 
speaking about has a strong dependence on the form of the Hamiltonian one is dealing 
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with. 

We proceed to randomly generate initial states |\E') as explained above. For each 
state we calculate: 

(i) its concurrence C^, and 

(ii) the first minimum that the fidelity attains during the time evolution of 
For a given pure state |$), tells us "how far" can |$) travel in 5, before 
starting backwards towards itself, as guided by the Hamiltonian. Deeper valleys 
may be reached later in the periodic time evolution, but we are interested only 
in the one that is reached first. 

(iii) the time r required for the state |\E') to reach the first fidelity minimum 

(iv) the time interval Bound that arises by the following process: fixing first an 
arbitrary F value, |\E') can evolve to states |$) such that the overlap (fT^ between 
|\E') and |$) attains this value. Some time interval to = VLh/e [cf. Eq. ()14|) ] is 
needed to reach each of these states. The minimum possible such interval is 
called and given by the bound (fTTjl . Notice that the fidelity value F 
reached in this time interval needs not correspond to any minimum of the fidelity. 

Thus, for each |\E') we compute its concurrence Cj,, the (first) minimum of the 
fidelity and the time r (in units of Ti Bound), i-e., the ratio t /T^ Bound- This 

allows us to build up an association connecting each |\E') to these three quantities: 

1^) ^ {C, F„„„, t/Tl. Bound}. (20) 

A representative group of those states for which -F^m G [0.35, 0.4] is depicted in Fig. 
1(b). Their behavior is quite different than those of Fig. 1(a). For these Fmin values 
there exist two different families of states, corresponding to the two extant possibilities 
for Tl Bound [cf. Eq. fITTj) ]. The lower one corresponds to those states for which the 
bound is determined by its expectation-energy value E [the first one in Eq. ()17p]. 
States for which the bound Bound is determined by its energy spread /S.E belong to 
the upper group. None of these states (for both sets) reach the bound r/Ti^Bound = 1- 
The maximum value for the ratio t/Tl, Bound of Fig. 1(a) is clearly exceeded here. 
There also exists a forbidden C zone for states with Fmin different than 0. To acquire 
a global perspective regarding these families of states, for any value of Fmin, "we will 
study the dependence of their fastest evolving states [Pi and P2 in Fig. 1(b)] with Fmin- 
These rapidly evolving states achieve the minimum of the fidelity in a time flmin = 7i/2. 
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0.2 0.4 0.6 0.8 1 

FIG. 2: Concurrence for the fastest pure states compatible with a given value of as 
given by Eq. (|22() . The horizontal and vertical lines cross at the point corresponding to the 
fastest states (points Pi and P2) of Fig. lb. 

Also, from Eq. ()15|). we realize that the minimum of the fidelity cannot be reached in 
a time shorted than Qmin = 7r/2. If proper account of the normalization (fTUj) of the 
initial state is taken, the only compatible state parametrization turns out to be 



|ci|' = |C2|' = 0, 
I 1 2 1 ~l~ vFrni 

Co = 7; 



2 1 — 



\cs\' = Y ■ ^21) 

Given such a parametrization and minding Eq. (jlll) . we ascertain that the con- 
currence Ci.Pi for these fast states is the same for our two families, being completely 
determined once the value of Fmin is fixed, 

CL.Pi = -x/l — Fmin- (22) 

This dependence on Fmin is illustrated if Fig. 2. For a given value of Fmin, the 
concurrence Cl p- of these fast states also coincides with the maximal concurrence 
value allowed for. Thus, determining the concurrence of the fastest states compatible 
with such fidelity is tantamount to finding the forbidden C zone for such Fmin value. 
Only in the special instance of states capable to reach (in their evolution-trajectory) 
an orthogonal counterpart (i.e., Fmin = 0) can we obtain any possible concurrence 
value. For the Fmin range depicted in Fig. 1(b), the points Pi and P2 correspond to 
Fmin = 0.35. According to Eq. their concurrence is C^.p. = 0.806, as shown in 
Figs. 1(b) and 2. 
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FIG. 3: r/r^.p. for the fastest pure states compatible with a given value of Fmin as given by 
Eq. PH) for t/Tl^p-^ and by Eq. (^5]) for t/Tl,p^. The horizontal and vertical lines cross at 
two different points corresponding to the fastest states (points Pi and P2) of Fig. 1(b). The 
upper crossing corresponds to P2 and the lower one to Pi. 

Once we know the concurrence value for these special states, we ascertain the time 
(in Tl Bound units) required to reach F^in- We need first to compute the relation 
between Fmin and either the mean energy E or the energy spread AE. Using the 
parametrization (pT| one obtains 



E = e{l-^/F~), 
AE = eVl - Fmin. (23) 

Since we know that the shortest possible time needed to reach this Fmin value is 
^min = 7r/2, using (fTTjl we easily find the ratio t/Tl, Bound corresponding to the "fastest 
states" of our two families. 



T 1 VFmin^ {24) 



mm 



Tl.P2 l3{Fmin) ^ ^ 

Both quantities are depicted in Fig. 3. The higher the Fmin value, the more apart 
the two families get. They only coincide in the Fmin = case, that is, for those 
states that evolve to an orthogonal one. We can apply these results to the special case 
considered in Fig. 1(b) {Fmin £ [0.35,0.4]), remembering that in such Fmin range the 
points Pi and P2 correspond to Fmin = 0.35. For the fast state corresponding to point 
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Pi we have E = 0.408 (in e units) and t/T^^p^ = 1.091. For the state corresponding 
to the point P2 one finds AE = 0.806 (in e units) and t/Tl,p^ = 1.351. 

By recourse to numerical simulation we have also found that the number of states 
that evolve according to Eqs. or (^3)) is a function of the value Emm- For = 
the bound for all the involved states is given by T^, Bound = P{E){nh/2AE) j2o|. For 
greater Emm values the situation changes. If this value is large enough, approximately 
half of the states belong to one of the families, while the rest are accrued to the 
remaining one. 



III. QUANTUM SPEED LIMIT FOR MIXED STATES 

Given an initial mixed state p{0) and using the Hamiltonian of the last section, we 
can easily calculate the corresponding density matrix of the system at a given time t 



Pit) 



( Pii Pi2e'^^^ Pne'^^^ pue^^^^ \ 

P3ie-*'^^ P32e'^~'' P33 P34e*'^^ 



(26) 



\ P4ie p42e ^ p43e " P44 / 

where pij = Pij{0). 

For mixed states the fidelity's expression adopts the well known expression 



F{p{0),pit)) = {Tr[^vW)pit)vW)]r- (27) 

In the case of pure states, this fidelity reduces to the probability (fT2j) . For the case 
treated here, determined by the Hamiltonian Hj, the fidelity for pure states is given 
by Eq. p5|) from which we realize that it oscillates in time. For mixed states such 
kind of analytical expression for the fidelity in not available, but one can compute the 
fidelity numerically and observe a similar behavior. 

To study the case of mixed states we follow the same methodology used in the 
previous section for pure states. Thus, we randomly generate states p in the two- 
qubits space of mixed states S (of 15 dimensions). We can thus classify the values of 
the ratio t/T^, Bound according to their corresponding values of the concurrence and 
the fidelity in such a mapping. We also fix our attention on the "concurrence fidelity" - 
plane. Our numerical-sampling procedure will start by constructing a fine enough grid 
in the (F, C) plane. We will have thus divided the plane into a large but discrete 
number of "windows". Each window, of course, contains many states p. We will assign 
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to all of them the same pair of {F, C) values. Notice that, for these distinct states p, 
r, and Tl, Bound will in general be different. We thus average over them, but omit, for 
notational simplicity's sake, the "< >" signs. The end-result is that we get a list of 
three quantities for each grid, namely, 

1. C, 

2. F, 

. Bound- 

We also store information regarding the times t at which the fidelity achieves some 
arbitrary fixed value, not necessarily connected with minima in any sense of the word. 
Specifically, for each p, we have selected intermediate values of the fidelity between its 
initial, and maximum, value F = 1, and its (first-) minimum value -F^j^, according to 
intervals of size 0.05 (F = 1, 0.95, 0.9, 0.85, . . .). For these fidelities, we have stored the 
associated quantities C, t/Tl, Bound- 
As stated, the time evolution for mixed states is of a periodic nature and the oscilla- 
tion strongly depends on the Hamiltonian form. One can circumvent to a considerable 
extent this F dependence on the form by using only high-fidelity values, for which, 
obviously, F minima cannot be reached in arbitrarily small time intervals. A major- 
ity of the states p G 5 attain these high-fidelity values but do not achieve, instead, 
lower ones. Using high-fidelity values is then tantamount to considering most of our 
randomly generated states G S. 

Let us focus our attention upon t/Tl, Bound- As mentioned before, the quantity 
Tl. Bound [cf. Eq. flTTjl ] (also known as the quantum speed hmit time) is the lower 
bound for the temporal interval required so as to evolve, from a state p to a state 
a, in such a manner that the pair (p, cr) of companion states corresponds to a given 
fidelity F. We will first study this quantity for fixed fidelity values. In Fig. 4 we 
plot the value of t/Tl, Bound vs the concurrence C for mixed states evolving to high 
fidelity companion-states. There exists a clear correlation between the quantum speed 
evolution time and the concurrence. The more entangled a state is, the less time it 
takes to reach a companion state such that the generalized overlap between them is 
F. Indeed, this time is seen to be close to the limit Tl. Bound for high C values. The 
relation between t/Tl. Bound and the concurrence C does not seem to strongly depend 
on the specific Fmm value considered. Contrariwise, for pure states it is only when 
we consider small fidelity values that a clear correlation between t/T^ Bound and the 
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FIG. 4: t/Tl. Bound for mixed states that evolve to several fixed {^F — 0.95 and F — 0.8) and 
minimum {Fmin = 0.925 and Fmin = 0.825) values of the fidelity. See text for details. 



concurrence C is observed. These small fixed values of the fidelity are obviously very 
close to its corresponding Fmin value. Remind that we have shown in the previous 
section that the correlation between t/Tl Bound and the concurrence C does exist for 
Fmin- If one selects pure states with high fixed values of the fidelity F, the correlation 
we are here speaking about is not detected, as t/Tl, Bound is approximately constant 
and close to unity for all possible values of the concurrence C. 

Thus far we have considered arbitrary fidelities, not necessarily linked to minima 
of this quantity. In Fig. 4 we also consider time intervals r needed to reach the first 
fidelity minimum. The resulting situation resembles the one just described above. The 
number of highly entangled (C > 0.9) mixed states evolving to companion states with 
a high value of the fidelity {F > 0.75) is very small. This entails that we cannot 
numerically obtain enough states in this zone to perform our averaging procedure in 
a reliable manner. Although the limit value t/Tl. Bound = 1 is not reached for maxi- 
mally entangled states, mixed states exhibit also in this case a nitid speed-concurrence 
correlation. Consequently, their evolution speed strongly depends on entanglement 
degree. 



IV. MAXIMALLY ENTANGLED MIXED STATES (MEMS) AND IH 
STATES 

As we are interested in putative relations between entanglement and the speeding 
up of the quantum evolution, the study of some special types of states should be of 
interest. In particular, the so-called MEMS l2l( are states that have the maximum 



14 



possible amount of entanglement of formation. We remark on the fact that MEM 
states have recently been encountered in the laboratory The associated 

density matrix is written in terms of a variable x that ranges in [0, 1]. In the basis 
referred to in the Sec. II their representative matrices read 



Pmems 



{2i 



( g{x) x/2\ 
l-2g{x) 


\x/2 g{x) ) 
with g{x) = 1/3 for < a; < 2/3, and g{x) = x/2 for 2/3<x < 1. 

Also of great interest are the so-called Ishizaka and Hiroshima (IH) states [2^, 
whose entanglement degree cannot be increased by acting on them with logic gates. 
Of course, MEMS are a special instance of the IH class. The associated pin matrices. 



of eigenvalues pi 



1,2,3,4), read 



PlH 



P2 



\ 



g P3+P1 P3-P1 Q 

£^ 2^ 



(29) 



V P4 y 

where the eigenvalues are size ordered: Pi > P2 > > P4- If one compares Eqs. ()29|) 
and ()26|) it is easy to see that IH states can "evolve" only if 6a ^ Sb- This entails that 
we cannot use here the same Hamiltonian (Hj) employed in the preceding sections (see 
Sec. II). For IH states we will use the values 6b = I and 6a = 2, namely, we employ 
a local Hamiltonian Hjj, whose diagonal is 0,e,2e,3e. For MEM states we will also 
use this Hamiltonian because we want to compare its associated results with those of 
the IH states. If one uses the hamiltonian Hj with MEM states the ensuing results 
resemble those of the preceding section. Thus, comparison can be made with the 
results of mixed states detailed there. In the MEMS instance the fidelity can readily 
be computed and reads 



with 



MEMS 




(30) 



A 



2g{x) 
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B = 45((x)^ + a;^cos(A+a), 

C = x\cos{A+a) -l)(8g{x)^ + x\cos{A+a) -l)y (31) 

The MEMS's fidelity expression (1311) also applies, with different coefficients, to the 
IH states case. Their corresponding coefficients are 



AlH=P2 +P4, 

BiH = {Pi + P^f + {Pi - P'if COS (A~a), 

Cm = (l + cos (A-«)) (^(pi^ + P32) {l + cos (A-^)) 

+ 2piP3(3-cos(A-a))^. (32) 

The oscillating part of the MEMS fidelity seems to depend on cos {5a + 5b)^ and 
the minimum fidelity value coincides with the minimum of its oscillating part, i.e., the 
minimum fidelity is achieved at Q^^^'^^ = vr/ (Sa + Sb)- For the IH-states the situation 
is similar, the fidelity depends on cos {6a — Sb)^ and its first minimum is achieved at 

We can also obtain an analytic expression for the expectation value of the hamilto- 
nian E and its fluctuation AE in the case of the MEMS states. We find 



Emems = e {S B + g{x)A ), 

AEmems = e \lg{x){6B^ + 5a') - g{xf{A-f. (33) 



For the Hamiltonian considered here E is always greater than AE. According to 
Eq. (|T7jl (see also the paragraph following it) one looks for the maximum of a pair 
of quantities. Here Tl, Bound is always equal to that one depending on AE, because 
(3{F) is always greater than or equal to a{F). Taking into account all these results, 
the equation for r turns out to read 

T ^'irii'^i AE (34) 



Tl. Bound e arccos ^/F ' 
and one must substitute AE and F with their pertinent associated values, depending 

on (i) which zone we are working in and (ii) which Hamiltonian we are referring to. 

For VLmin we obtain the values 7r/2, for {5b = 1 , = 1), and 7r/3, for {5b = 2 , 

5a = 1). 
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FIG. 5: t/Tl. Bound for MEMS with hamiltonian Hn and the corresponding average for the 
IH case. See text for details 

For the IH states we obtain the following equations for E and A£' 

EiH = |(piA- + (p3 + 2p4)A+), 
^EiH = |^2piA+A- + 2(p3 + 2p4)(A+)' - (piA" + (pa + 2p4)A+)'. (35) 

Notice that the MEM states are completely determined by the parameter which 
corresponds to the value of the concurrence {C = x). For a given value of the concur- 
rence C there exists only one value for that magnitudes we are interested on: Fmim 
E, AE and TL.Bound- This means that if we want to analyze these states we must do 
it using F^in instead of fixed, arbitrary values of the fidelity. If we use fixed, arbi- 
trary fidelity values K, one can always detect a range of C values for which no MEMS 
characterized by K exist. This feature constitutes a great difference with respect to 
what happens for the general mixed states case discussed above, where one has many 
different states compatible with a given value of the concurrence C. Such is the case 
for IH states, for which we can average such magnitudes as we did earlier for general 
mixed states in Sec. III. We must mind this difference between MEMS and IH states 
if we want to compare MEM results to IH ones. As stated before, for an arbitrary 
MEM state the concurrence, say C = K is fixed, and so is the value of F„iin. But for 
the same K value there exist many IH states characterized by a wide range of possible 
Fmin values. 

In order to be able to compare IH and MEM states we have used the following 
criterion: consider those IH states of concurrence C — K with an value greater 
than that pertaining to the associated K MEM state. In Fig. 5 we depict t/Tl, Bound 
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vs C for MEMS and for those IH states that fulfill the above criterion. The behavior of 
separable states, with exceedingly large t/Tl. Bound values, is very different from that 
of highly entangled ones. The ensuing differences are larger than for the general mixed 
states studied in the preceding section. There is a clear difference between the two 
MEM zones arising out of the x— MEM parametrization. For weakly entangled states, 
't/Tl. Bound achieves very high values. For highly entangled states (the second MEM 
zone), the situation is the opposite. In this last zone t/Tl. Bound tends to saturate its 
lower bound. The IH states have a similar behavior than the MEMS, although its 
corresponding average of the ratio t/Tl, Bound is always greater than the corresponding 
MEMS ratio. 

V. CONCLUSIONS 

For mixed states of bipartite systems ruled by a general local Hamiltonian we have 
put forward rather convincing evidence of the clear correlation extant between concur- 
rence and speed of time evolution. The more entangled an initial state pi is, the less 
time (in units of Tl. Bound) it takes to evolve to another state p2, no matter what the pi 
-p2 fidehty is. In the case of pure states the correlation is strong for states that evolve 
to a minimum of the fidelity. For some special mixed states, namely, the so-called 
maximally entangled mixed states, the correlation between concurrence and the speed 
of time evolution becomes more acute than in the case of general mixed states. 
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